Local and systemic replicative fitness for viruses in specialist, generalist, and non-specialist interactions with salmonid hosts.

Host tissues represent diverse resources or barriers for pathogen replicative fitness. We tested whether viruses in specialist, generalist, and non-specialist interactions replicate differently in local entry tissue (fin), and systemic target tissue (kidney) using infectious hematopoietic necrosis virus (IHNV) and three salmonid fish hosts. Virus tissue replication was host specific, but one feature was shared by specialists and the generalist which was uncommon in the non-specialist interactions: high host entry and replication capacity in the local tissue after contact. Moreover, specialists showed increased replication in systemic target tissues early after host contact. By comparing ancestral and derived IHNV viruses, we also characterized replication tradeoffs associated with specialist and generalist evolution. Compared with the ancestral virus, a derived specialist gained early local replicative fitness in the new host but lost replicative fitness in the ancestral host. By contrast, a derived generalist showed small replication losses relative to the ancestral virus in the ancestral host but increased early replication in the local tissue of novel hosts. This study shows that the mechanisms of specialism and generalism are host specific and that local and systemic replication can contribute differently to overall within host replicative fitness for specialist and generalist viruses.

Assortative mating for reproductive timing affects population recruitment and resilience in a quantitative genetic model.

Quantitative models that simulate the inheritance and evolution of fitness-linked traits offer a method for predicting how environmental or anthropogenic perturbations can affect the dynamics of wild populations. Random mating between individuals within populations is a key assumption of many such models used in conservation and management to predict the impacts of proposed management or conservation actions. However, recent evidence suggests that non-random mating may be underestimated in wild populations and play an important role in diversity-stability relationships. Here we introduce a novel individual-based quantitative genetic model that incorporates assortative mating for reproductive timing, a defining attribute of many aggregate breeding species. We demonstrate the utility of this framework by simulating a generalized salmonid lifecycle, varying input parameters, and comparing model outputs to theoretical expectations for several eco-evolutionary, population dynamic scenarios. Simulations with assortative mating systems resulted in more resilient and productive populations than those that were randomly mating. In accordance with established ecological and evolutionary theory, we also found that decreasing the magnitude of trait correlations, environmental variability, and strength of selection each had a positive effect on population growth. Our model is constructed in a modular framework so that future components can be easily added to address pressing issues such as the effects of supportive breeding, variable age structure, differential selection by sex or age, and fishery interactions on population growth and resilience. With code published in a public Github repository, model outputs may easily be tailored to specific study systems by parameterizing with empirically generated values from long-term ecological monitoring programs.

Population-size history inferences from the coho salmon (Oncorhynchus kisutch) genome.

Coho salmon (Oncorhynchus kisutch) are a culturally and economically important species that return from multiyear ocean migrations to spawn in rivers that flow to the Northern Pacific Ocean. Southern stocks of coho salmon in Canada and the United States have significantly declined over the past quarter century, and unfortunately, conservation efforts have not reversed this trend. To assist in stock management and conservation efforts, we generated a chromosome-level genome assembly. We also resequenced the genomes of 83 coho salmon across the North American range to identify nucleotide variants and understand the demographic histories of these salmon by modeling effective population size from genome-wide data. From demographic history modeling, we observed reductions in effective population sizes between 3,750 and 8,000 years ago for several northern sampling sites, which may correspond to bottleneck events during recolonization after glacial retreat.

Long-distance migration is a major factor driving local adaptation at continental scale in Coho salmon.

Inferring the genomic basis of local adaptation is a long-standing goal of evolutionary biology. Beyond its fundamental evolutionary implications, such knowledge can guide conservation decisions for populations of conservation and management concern. Here, we investigated the genomic basis of local adaptation in the Coho salmon (Oncorhynchus kisutch) across its entire North American range. We hypothesized that extensive spatial variation in environmental conditions and the species' homing behaviour may promote the establishment of local adaptation. We genotyped 7829 individuals representing 217 sampling locations at more than 100,000 high-quality RADseq loci to investigate how recombination might affect the detection of loci putatively under selection and took advantage of the precise description of the demographic history of the species from our previous work to draw accurate population genomic inferences about local adaptation. The results indicated that genetic differentiation scans and genetic-environment association analyses were both significantly affected by variation in recombination rate as low recombination regions displayed an increased number of outliers. By taking these confounding factors into consideration, we revealed that migration distance was the primary selective factor driving local adaptation and partial parallel divergence among distant populations. Moreover, we identified several candidate single nucleotide polymorphisms associated with long-distance migration and altitude including a gene known to be involved in adaptation to altitude in other species. The evolutionary implications of our findings are discussed along with conservation applications.

Variation in within-host replication kinetics among virus genotypes provides evidence of specialist and generalist infection strategies across three salmonid host species.

Theory of the evolution of pathogen specialization suggests that a specialist pathogen gains high fitness in one host, but this comes with fitness loss in other hosts. By contrast, a generalist pathogen does not achieve high fitness in any host, but gains ecological fitness by exploiting different hosts, and has higher fitness than specialists in nonspecialized hosts. As a result, specialist pathogens are predicted to have greater variation in fitness across hosts, and generalists would have lower fitness variation across hosts. We test these hypotheses by measuring pathogen replicative fitness as within-host viral loads from the onset of infection to the beginning of virus clearance, using the rhabdovirus infectious hematopoietic necrosis virus (IHNV) in salmonid fish. Based on field prevalence and virulence studies, the IHNV subgroups UP, MD, and L are specialists, causing infection and mortality in sockeye salmon, steelhead, and Chinook salmon juveniles, respectively. The UC subgroup evolved naturally from a UP ancestor and is a generalist infecting all three host species but without causing severe disease. We show that the specialist subgroups had the highest peak and mean viral loads in the hosts in which they are specialized, and they had low viral loads in nonspecialized hosts, resulting in large variation in viral load across hosts. Viral kinetics show that the mechanisms of specialization involve the ability to both maximize early virus replication and avoid clearance at later times, with different mechanisms of specialization evident in different host-virus combinations. Additional nuances in the data included different fitness levels for nonspecialist interactions, reflecting different trade-offs for specialist viruses in other hosts. The generalist UC subgroup reached intermediate viral loads in all hosts and showed the smallest variation in fitness across hosts. The evolution of the UC generalist from an ancestral UP sockeye specialist was associated with fitness increases in steelhead and Chinook salmon, but only slight decreases in fitness in sockeye salmon, consistent with low- or no-cost generalism. Our results support major elements of the specialist-generalist theory, providing evidence of a specialist-generalist continuum in a vertebrate pathogen. These results also quantify within-host replicative fitness trade-offs resulting from the natural evolution of specialist and generalist virus lineages in multi-host ecosystems.

Temperature Variation and Host Immunity Regulate Viral Persistence in a Salmonid Host.

Environmental variation has important effects on host-pathogen interactions, affecting large-scale ecological processes such as the severity and frequency of epidemics. However, less is known about how the environment interacts with host immunity to modulate virus fitness within hosts. Here, we studied the interaction between host immune responses and water temperature on the long-term persistence of a model vertebrate virus, infectious hematopoietic necrosis virus (IHNV) in steelhead trout (Oncorhynchus mykiss). We first used cell culture methods to factor out strong host immune responses, allowing us to test the effect of temperature on viral replication. We found that 15 ∘C water temperature accelerated IHNV replication compared to the colder 10 and 8 ∘C temperatures. We then conducted in vivo experiments to quantify the effect of 6, 10, and 15 ∘C water temperatures on IHNV persistence over 8 months. Fish held at 15 and 10 ∘C were found to have higher prevalence of neutralizing antibodies compared to fish held at 6 ∘C. We found that IHNV persisted for a shorter time at warmer temperatures and resulted in an overall lower fish mortality compared to colder temperatures. These results support the hypothesis that temperature and host immune responses interact to modulate virus persistence within hosts. When immune responses were minimized (i.e., in vitro) virus replication was higher at warmer temperatures. However, with a full potential for host immune responses (i.e., in vivo experiments) longer virus persistence and higher long-term virulence was favored in colder temperatures. We also found that the viral RNA that persisted at later time points (179 and 270 days post-exposure) was mostly localized in the kidney and spleen tissues. These tissues are composed of hematopoietic cells that are favored targets of the virus. By partitioning the effect of temperature on host and pathogen responses, our results help to better understand environmental drivers of host-pathogen interactions within hosts, providing insights into potential host-pathogen responses to climate change.

Virulence and Infectivity of UC, MD, and L Strains of Infectious Hematopoietic Necrosis Virus (IHNV) in Four Populations of Columbia River Basin Chinook Salmon.

Infectious Hematopoietic Necrosis Virus (IHNV) infects juvenile salmonid fish in conservation hatcheries and aquaculture facilities, and in some cases, causes lethal disease. This study assesses intra-specific variation in the IHNV susceptibility of Chinook salmon (Oncorhynchus tshawytscha) in the Columbia River Basin (CRB), in the northwestern United States. The virulence and infectivity of IHNV strains from three divergent virus genogroups are measured in four Chinook salmon populations, including spring-run and fall-run fish from the lower or upper regions of the CRB. Following controlled laboratory exposures, our results show that the positive control L strain had significantly higher virulence, and the UC and MD strains that predominate in the CRB had equivalently low virulence, consistent with field observations. By several experimental measures, there was little variation in host susceptibility to infection or disease. However, a small number of exceptions suggested that the lower CRB spring-run Chinook salmon population may be less susceptible than other populations tested. The UC and MD viruses did not differ in infectivity, indicating that the observed asymmetric field prevalence in which IHNV detected in CRB Chinook salmon is 83% UC and 17% MD is not due to the UC virus being more infectious. Overall, we report little intra-species variation in CRB Chinook salmon susceptibility to UC or MD IHNV infection or disease, and suggest that other factors may instead influence the ecology of IHNV in the CRB.

Correction: Demographic history shaped geographical patterns of deleterious mutation load in a broadly distributed Pacific Salmon.

[This corrects the article DOI: 10.1371/journal.pgen.1008348.].

Heterogeneous genetic basis of age at maturity in salmonid fishes.

Understanding the genetic basis of repeated evolution of the same phenotype across taxa is a fundamental aim in evolutionary biology and has applications in conservation and management. However, the extent to which interspecific life-history trait polymorphisms share evolutionary pathways remains underexplored. Here, we address this gap by studying the genetic basis of a key life-history trait, age at maturity, in four species of Pacific salmonids (genus Oncorhynchus) that exhibit intra- and interspecific variation in this trait-Chinook Salmon, Coho Salmon, Sockeye Salmon, and Steelhead Trout. We tested for associations in all four species between age at maturity and two genome regions, six6 and vgll3, that are strongly associated with the same trait in Atlantic Salmon (Salmo salar). We also conducted a genome-wide association analysis in Steelhead to assess whether additional regions were associated with this trait. We found the genetic basis of age at maturity to be heterogeneous across salmonid species. Significant associations between six6 and age at maturity were observed in two of the four species, Sockeye and Steelhead, with the association in Steelhead being particularly strong in both sexes (p = 4.46 × 10-9 after adjusting for genomic inflation). However, no significant associations were detected between age at maturity and the vgll3 genome region in any of the species, despite its strong association with the same trait in Atlantic Salmon. We discuss possible explanations for the heterogeneous nature of the genetic architecture of this key life-history trait, as well as the implications of our findings for conservation and management.

Power of a dual-use SNP panel for pedigree reconstruction and population assignment.

The use of high-throughput, low-density sequencing approaches has dramatically increased in recent years in studies of eco-evolutionary processes in wild populations and domestication in commercial aquaculture. Most of these studies focus on identifying panels of SNP loci for a single downstream application, whereas there have been few studies examining the trade-offs for selecting panels of markers for use in multiple applications. Here, we detail the use of a bioinformatic workflow for the development of a dual-purpose SNP panel for parentage and population assignment, which included identifying putative SNP loci, filtering for the most informative loci for the two tasks, designing effective multiplex PCR primers, optimizing the SNP panel for performance, and performing quality control steps for downstream applications. We applied this workflow to two adjacent Alaskan Sockeye Salmon populations and identified a GTseq panel of 142 SNP loci for parentage and 35 SNP loci for population assignment. Only 50-75 panel loci were necessary for >95% accurate parentage, whereas population assignment success, with all 172 panel loci, ranged from 93.9% to 96.2%. Finally, we discuss the trade-offs and complexities of the decision-making process that drives SNP panel development, optimization, and testing.

Infectious hematopoietic necrosis virus specialization in a multihost salmonid system.

Many pathogens interact and evolve in communities where more than one host species is present, yet our understanding of host-pathogen specialization is mostly informed by laboratory studies with single species. Managing diseases in the wild, however, requires understanding how host-pathogen specialization affects hosts in diverse communities. Juvenile salmonid mortality in hatcheries caused by infectious hematopoietic necrosis virus (IHNV) has important implications for salmonid conservation programs. Here, we evaluate evidence for IHNV specialization on three salmonid hosts and assess how this influences intra- and interspecific transmission in hatchery-reared salmonids. We expect that while more generalist viral lineages should pose an equal risk of infection across host types, viral specialization will increase intraspecific transmission. We used Bayesian models and data from 24 hatcheries in the Columbia River Basin to reconstruct the exposure history of hatcheries with two IHNV lineages, MD and UC, allowing us to estimate the probability of juvenile infection with these lineages in three salmonid host types. Our results show that lineage MD is specialized on steelhead trout and perhaps rainbow trout (both Oncorhynchus mykiss), whereas lineage UC displayed a generalist phenotype across steelhead trout, rainbow trout, and Chinook salmon. Furthermore, our results suggest the presence of specialist-generalist trade-offs because, while lineage UC had moderate probabilities of infection across host types, lineage MD had a small probability of infection in its nonadapted host type, Chinook salmon. Thus, in addition to quantifying probabilities of infection of socially and economically important salmonid hosts with different IHNV lineages, our results provide insights into the trade-offs that viral lineages incur in multihost communities. Our results suggest that knowledge of the specialist/generalist strategies of circulating viral lineages could be useful in salmonid conservation programs to control disease.

Demographic history shaped geographical patterns of deleterious mutation load in a broadly distributed Pacific Salmon.

A thorough reconstruction of historical processes is essential for a comprehensive understanding of the mechanisms shaping patterns of genetic diversity. Indeed, past and current conditions influencing effective population size have important evolutionary implications for the efficacy of selection, increased accumulation of deleterious mutations, and loss of adaptive potential. Here, we gather extensive genome-wide data that represent the extant diversity of the Coho salmon (Oncorhynchus kisutch) to address two objectives. We demonstrate that a single glacial refugium is the source of most of the present-day genetic diversity, with detectable inputs from a putative secondary micro-refugium. We found statistical support for a scenario whereby ancestral populations located south of the ice sheets expanded recently, swamping out most of the diversity from other putative micro-refugia. Demographic inferences revealed that genetic diversity was also affected by linked selection in large parts of the genome. Moreover, we demonstrate that the recent demographic history of this species generated regional differences in the load of deleterious mutations among populations, a finding that mirrors recent results from human populations and provides increased support for models of expansion load. We propose that insights from these historical inferences should be better integrated in conservation planning of wild organisms, which currently focuses largely on neutral genetic diversity and local adaptation, with the role of potentially maladaptive variation being generally ignored.

Genomic and phenotypic effects of inbreeding across two different hatchery management regimes in Chinook salmon.

Genomic approaches permit direct estimation of inbreeding and its effect on fitness. We used genomic-based estimates of inbreeding to investigate their relationship with eight adult traits in a captive-reared Pacific salmonid that is released into the wild. Estimates were also used to determine whether alternative broodstock management approaches reduced risks of inbreeding. Specifically, 1,100 unlinked restriction-site associated (RAD) loci were used to compare pairwise relatedness, derived from a relationship matrix, and individual inbreeding, estimated by comparing observed and expected homozygosity, across four generations in two hatchery lines of Chinook salmon that were derived from the same source. The lines are managed as "integrated" with the founding wild stock, with ongoing gene flow, and as "segregated" with no gene flow. While relatedness and inbreeding increased in the first generation of both lines, possibly due to population subdivision caused by hatchery initiation, the integrated line had significantly lower levels in some subsequent generations (relatedness: F2 -F4 ; inbreeding F2 ). Generally, inbreeding was similar between the lines despite large differences in effective numbers of breeders. Inbreeding did not affect fecundity, reproductive effort, return timing, fork length, weight, condition factor, and daily growth coefficient. However, it delayed spawn timing by 1.75 days per one standard deviation increase in F (~0.16). The results indicate that integrated management may reduce inbreeding but also suggest that it is relatively low in a small, segregated hatchery population that maximized number of breeders. Our findings demonstrate the utility of genomics to monitor inbreeding under alternative management strategies in captive breeding programs.

Publisher Correction: Sex-dependent dominance maintains migration supergene in rainbow trout.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Conservation and Management of Salmon in the Age of Genomics.

Salmon were among the first nonmodel species for which systematic population genetic studies of natural populations were conducted, often to support management and conservation. The genomics revolution has improved our understanding of the evolutionary ecology of salmon in two major ways: (a) Large increases in the numbers of genetic markers (from dozens to 104-106) provide greater power for traditional analyses, such as the delineation of population structure, hybridization, and population assignment, and (b) qualitatively new insights that were not possible with traditional genetic methods can be achieved by leveraging detailed information about the structure and function of the genome. Studies of the first type have been more common to date, largely because it has taken time for the necessary tools to be developed to fully understand the complex salmon genome. We expect that the next decade will witness many new studies that take full advantage of salmonid genomic resources.

Sex-dependent dominance maintains migration supergene in rainbow trout.

Males and females often differ in their fitness optima for shared traits that have a shared genetic basis, leading to sexual conflict. Morphologically differentiated sex chromosomes can resolve this conflict and protect sexually antagonistic variation, but they accumulate deleterious mutations. However, how sexual conflict is resolved in species that lack differentiated sex chromosomes is largely unknown. Here we present a chromosome-anchored genome assembly for rainbow trout (Oncorhynchus mykiss) and characterize a 55-Mb double-inversion supergene that mediates sex-specific migratory tendency through sex-dependent dominance reversal, an alternative mechanism for resolving sexual conflict. The double inversion contains key photosensory, circadian rhythm, adiposity and sex-related genes and displays a latitudinal frequency cline, indicating environmentally dependent selection. Our results show sex-dependent dominance reversal across a large autosomal supergene, a mechanism for sexual conflict resolution capable of protecting sexually antagonistic variation while avoiding the homozygous lethality and deleterious mutations associated with typical heteromorphic sex chromosomes.

Influence of a growth hormone transgene on the genetic architecture of growth-related traits: A comparative analysis between transgenic and wild-type coho salmon.

Genetic engineering has been increasingly applied to many commercially important plant and animal species, generating phenotypic changes that are not observed in natural populations and creating genetic interactions that have not experienced natural selection. The degree to and way in which such human-induced genetic variation interacts with the rest of the genome is currently largely unknown. Integrating such information into ecological and risk assessment frameworks is crucial to understand the potential effects of genetically modified organisms in natural environments. Here, we performed QTL mapping to investigate the genetic architecture of growth-related traits in nontransgenic (NT) and growth hormone transgenic (T) coho salmon with large changes in growth and related physiology, with the aim of identifying how an inserted transgene might influence the opportunity for selection. These fish shared the same parental genetic background, thus allowing us to determine whether the same or different loci influence these traits within the two groups. The use of over 1,700 loci, derived from restriction site-associated DNA sequencing, revealed that different genomic regions were linked with growth over time between the two groups. Additionally, the effect sizes of detected QTL appear to have been influenced by the transgene. Direct comparison of QTL between the T and NT fish during two size-matched periods identified little overlap in their location. Taken together, the results showed that the transgene altered the genetic basis of growth-related traits in this species. The study has important implications for effective conservation and management of wild populations experiencing introduction of transgenes. Evolutionary changes and their ecological consequences may occur at different rates and in different directions in NT versus T individuals in response to selection. Thus, assessments of phenotypic change, and hence ecological risk, should be determined periodically to evaluate whether initial estimates made with founder strains remain valid.

Genomewide association analyses of fitness traits in captive-reared Chinook salmon: Applications in evaluating conservation strategies.

A novel application of genomewide association analyses is to use trait-associated loci to monitor the effects of conservation strategies on potentially adaptive genetic variation. Comparisons of fitness between captive- and wild-origin individuals, for example, do not reveal how captive rearing affects genetic variation underlying fitness traits or which traits are most susceptible to domestication selection. Here, we used data collected across four generations to identify loci associated with six traits in adult Chinook salmon (Oncorhynchus tshawytscha) and then determined how two alternative management approaches for captive rearing affected variation at these loci. Loci associated with date of return to freshwater spawning grounds (return timing), length and weight at return, age at maturity, spawn timing, and daily growth coefficient were identified using 9108 restriction site-associated markers and random forest, an approach suitable for polygenic traits. Mapping of trait-associated loci, gene annotations, and integration of results across multiple studies revealed candidate regions involved in several fitness-related traits. Genotypes at trait-associated loci were then compared between two hatchery populations that were derived from the same source but are now managed as separate lines, one integrated with and one segregated from the wild population. While no broad-scale change was detected across four generations, there were numerous regions where trait-associated loci overlapped with signatures of adaptive divergence previously identified in the two lines. Many regions, primarily with loci linked to return and spawn timing, were either unique to or more divergent in the segregated line, suggesting that these traits may be responding to domestication selection. This study is one of the first to utilize genomic approaches to demonstrate the effectiveness of a conservation strategy, managed gene flow, on trait-associated-and potentially adaptive-loci. The results will promote the development of trait-specific tools to better monitor genetic change in captive and wild populations.

A New Single Nucleotide Polymorphism Database for Rainbow Trout Generated Through Whole Genome Resequencing.

Single-nucleotide polymorphisms (SNPs) are highly abundant markers, which are broadly distributed in animal genomes. For rainbow trout (Oncorhynchus mykiss), SNP discovery has been previously done through sequencing of restriction-site associated DNA (RAD) libraries, reduced representation libraries (RRL) and RNA sequencing. Recently we have performed high coverage whole genome resequencing with 61 unrelated samples, representing a wide range of rainbow trout and steelhead populations, with 49 new samples added to 12 aquaculture samples from AquaGen (Norway) that we previously used for SNP discovery. Of the 49 new samples, 11 were double-haploid lines from Washington State University (WSU) and 38 represented wild and hatchery populations from a wide range of geographic distribution and with divergent migratory phenotypes. We then mapped the sequences to the new rainbow trout reference genome assembly (GCA_002163495.1) which is based on the Swanson YY doubled haploid line. Variant calling was conducted with FreeBayes and SAMtools mpileup, followed by filtering of SNPs based on quality score, sequence complexity, read depth on the locus, and number of genotyped samples. Results from the two variant calling programs were compared and genotypes of the double haploid samples were used for detecting and filtering putative paralogous sequence variants (PSVs) and multi-sequence variants (MSVs). Overall, 30,302,087 SNPs were identified on the rainbow trout genome 29 chromosomes and 1,139,018 on unplaced scaffolds, with 4,042,723 SNPs having high minor allele frequency (MAF > 0.25). The average SNP density on the chromosomes was one SNP per 64 bp, or 15.6 SNPs per 1 kb. Results from the phylogenetic analysis that we conducted indicate that the SNP markers contain enough population-specific polymorphisms for recovering population relationships despite the small sample size used. Intra-Population polymorphism assessment revealed high level of polymorphism and heterozygosity within each population. We also provide functional annotation based on the genome position of each SNP and evaluate the use of clonal lines for filtering of PSVs and MSVs. These SNPs form a new database, which provides an important resource for a new high density SNP array design and for other SNP genotyping platforms used for genetic and genomics studies of this iconic salmonid fish species.

A practical introduction to Random Forest for genetic association studies in ecology and evolution.

Large genomic studies are becoming increasingly common with advances in sequencing technology, and our ability to understand how genomic variation influences phenotypic variation between individuals has never been greater. The exploration of such relationships first requires the identification of associations between molecular markers and phenotypes. Here, we explore the use of Random Forest (RF), a powerful machine-learning algorithm, in genomic studies to discern loci underlying both discrete and quantitative traits, particularly when studying wild or nonmodel organisms. RF is becoming increasingly used in ecological and population genetics because, unlike traditional methods, it can efficiently analyse thousands of loci simultaneously and account for nonadditive interactions. However, understanding both the power and limitations of Random Forest is important for its proper implementation and the interpretation of results. We therefore provide a practical introduction to the algorithm and its use for identifying associations between molecular markers and phenotypes, discussing such topics as data limitations, algorithm initiation and optimization, as well as interpretation. We also provide short R tutorials as examples, with the aim of providing a guide to the implementation of the algorithm. Topics discussed here are intended to serve as an entry point for molecular ecologists interested in employing Random Forest to identify trait associations in genomic data sets.